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ABSTRACT 

An  analogy  between  damage  inflicted  by  a  single  point  hard  particle 
impact  and  damage  inflicted  by  inserting  a  flaw  of  known  dimensions  in  a 
static  tensile  coupon  is  discussed.  The  results  suggest  that  residual 
strength  can  be  predicted  as  a  function  of  kinetic  energy  of  impact  by  exe¬ 
cuting  two  experiments,  a  static  tensile  test  on  an  unflawed  specimen  and  a 
static  tensile  test  on  a  coupon  previously  subjected  to  a  single  point  impact. 
The  model  appears  to  be  accurate  for  impact  velocities  which  are  less  than 
the  penetration  velocity.  For  velocities  above  complete  penetration,  the 
residual  strength  is  identical  to  the  static  strength  of  a  coupon  with  a  hole 
having  the  same  diameter  as  the  impacting  particle.  Comparison  of  various 
materials  indicates  that  the  impact  strength  of  composite  materials  is 
strongly  influenced  by  the  strain  energy  to  failure  of  the  reinforcement. 


Accession  "or 

ins  i 

irxc  ?••; 

Unaanc  :  :  -ad  □ 

Justification _ 


By - 

Distribution/ 

Availability  Codes 
Av-  .1  and/or 
Dlst  .  (ji-olr.i 


iii 


«... 


AFML-TR-73-309 


TABLE  OF  CONTENTS 


SECTION 

I  INTRODUCTION 

II  ANALYSIS 

III  EXPERIMENTAL  PROCEDURE 

IV  DISCUSSION 

V  CONCLUSIONS 

APPENDIX  -  STIFFNESS  AND  FRACTURE  PREDICTIONS 
USING  LAMINATION  THEORY  AND  MAXIMUM 
STRAIN  FAILURE  CRITERION 


REFERENCES 


AFML-TR-7 3-309 


ILLUSTRATIONS 


FIGURE  PAGE 

1.  Comparison  Between  Strength  Reduction  Produced 

By  A  Drilled  Hole  and  the  Strength  Reduction 
Produced  By  A  Bullet  31 

2.  Residual  Strength  As  A  Function  of  Kinetic  Energy  32 

3.  Schematic  of  Residual  Strength  Analogy  33 

4.  Schematic  of  Experimental  Set-up  34 

5.  Effect  of  Matrix  on  Impact  Residual  Strength 

(MOD  II  Fiber)  35 

6.  Effect  of  Fiber  on  Impact  Residual  Strength 
(Fiber  Data  in  Table  1,  ERL  4617  Resin  System 

Used  in  All  Cases)  36 

7.  Effect  of  Specimen  Size  on  Impact  Residual  Strength 

(A-S/4617  Composites)  37 

8.  Effect  of  Specimen  Width  to  Projectile  Diameter 

on  Impact  Residual  Strength  (Scotchply/1002 
Composites)  38 

9.  Effect  of  Specimen  Width  to  Projectile  Diameter 

on  the  Value  of  K  (Scotchply/1002  Composites)  39 


vi 


i 


I 


AFML-TR-73-309 


SECTION  I 


INTRODUCTION 


Static  residual  strength  is  usually  defined  as  the  failure  stress  in  a 
uniaxial  tension  test  performed  after  the  specimen  has  been  subjected  to 
some  previous  load  history.  Such  information  is  a  useful  measure  of  dam¬ 
age  and  can  be  used  in  predicting  the  life  expectancy  of  a  composite  material 
subjected  to  a  specified  load  history  [l].  In  a  similar  manner,  residual 
strength  should  be  a  useful  measure  of  damage  in  a  composite  subjected  to 
a  local  impact  load,  as  well  as  a  means  of  screening  materials  for  potential 
application  in  structures  subjected  to  an  impact  load  such  as  turbine  engine 
fan  blades. 

It  has  been  previously  noted  [2]  that  a  hole  in  a  composite  material 
inflicted  by  a  bullet  produces  the  same  reduction  in  static  strength  as  a 
drilled  hole  of  the  same  diameter  as  the  bullet.  This  is  illustrated  in 
Figure  1  for  boron-epoxy,  graphite -epoxy,  and  glass -epoxy  laminates. 

Thus,  an  example  exists  in  which  the  damage  inflicted  by  a  local  hard  par¬ 
ticle  impact  can  be  equated,  in  terms  of  residual  strength,  to  damage 
inflicted  by  an  artificially  implanted  stress  concentration. 

Now  consider  a  tensile  coupon  subjected  to  a  localized  hard  particle 
impact  at  a  velocity  which  is  less  than  the  penetration  velocity.  If  this 
specimen  is  then  loaded  to  failure,  the  resulting  strength  will  be  less  than 
the  original  tensile  strength.  This  is  exactly  the  same  result 
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produced  by  implanting  a  small  through- the -thickness  crack  in  a  tensile 
coupon  and  then  stressing  it  to  failure.  Thus,  if  a  procedure  can  be 
developed  for  converting  the  impact  damage,  in  terms  of  residual  strength, 
to  an  equivalent  crack  of  known  dimensions,  the  residual  strength  can  be 
analyzed  in  terms  of  a  current  fracture  mechanics  model.  Such  an  analogy 
is  developed  and  executed  in  the  present  paper.  Comparison  between  the 
model  and  experimental  data  shows  good  agreement.  The  model  is  also 
used  to  compare  the  local  impact  resistance  of  various  types  of  composite 
materials . 
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SECTION  II 
ANALYSIS 

Consider  a  plate  of  orthotropic  construction  containing  a  narrow  slit 
of  length  2c  perpendicular  to  a  uniform  tensile  stress,  a,  applied  at 

infinity.  The  critical  strain  energy  release  rate,  Gjc,  has  been  deter- 

3 

mined  by  Sih,  Paris,  and  Irwin  [3] 


where  Kjc  is  the  critical  stress  intensity  factor  and  Sjj  are  orthotropic 
plate  compliances.  For  a  slit  in  an  orthotropic  plate,  the  value  of  Kic 
is  the  same  as  for  an  isotropic  material  [3],  i.e., 

Klc  =  (2) 


Substituting  Eq.  (2)  into  Eq.  (1)  yields  the  result 


For  an  isotropic  material  Eq.  (4)  reduces  to 

A  =  2ir  (5) 

3 

In  Reference  3,  the  Xj  axis  is  parallel  to  the  crack. 
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A  cursory  examination  of  Eq.  (3)  reveals  that  the  last  two  terms  on  the 
right-hand  side  constitute  the  work/unit  volume,  Wb,  necessary  to  break 
the  specimen,  i.e.,  the  area  under  the  stress-strain  curve  as  determined 
from  the  strain  field  at  a  distance  from  the  slit.  Thus, 

Glc  =  AcWb  (6) 

For  an  isotropic  material,  Eq.  (5)  in  conjunction  with  Eq.  (6)  yields 

Glc  =  2ircWb  (7) 

It  is  interesting  to  note  that  Rivlin  and  Thomas  [4]  have  shown  that  Eq. 

(7)  is  applicable  to  the  small  deformations  of  isotropic  materials  with 
nonlinear  constitutive  relations,  the  nonlinearity  being  accounted  for  by 
the  strain  energy  density  term.  Thus,  it  is  strongly  suspected  that  Eq.  (6) 
is  applicable  to  orthotropic  materials  with  nonlinear  constitutive  relations. 
This  means  that  Eq.  (6)  represents  a  very  useful  form  of  the  strain  energy 
release  rate. 

Now  let  us  assume  that  a  damage  zone  adjacent  to  a  stress  concentration 
constitutes  a  characteristic  volume  of  material  which  must  be  stressed  to 
a  critical  level  before  fracture.  Physically  such  a  zone  represents  an 
area  of  crazing  and  delamination  in  a  composite  material,  as  discussed 
previously  by  Halpin,  Jerina,  and  Johnson  [1].  Furthermore,  it  is 
assumed  that  this  volume  of  material  is  identified  by  a  characteristic 
dimension,  or  effective  flaw,  cQ.  In  the  absence  of  a  mechanically 
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implanted  flaw,  it  is  assumed  that  this  effective  flaw  determines  the 
strength  of  laminated  composites.  Under  these  assumptions  Eq.  (6) 
can  also  be  written  in  the  form 


G 


lc 


Ac  W 
o 


s 


(8) 


when  Wg  is  the  energy  under  the  stress- strain  curve  for  statically 
loaded  composites  without  a  mechanically  implanted  flaw. 

Using  the  damage  zone,  cQ,  in  the  same  manner  as  the  plastic  zone 
approximation  of  Irwin  [5],  an  effective  half  crack  length  c  +  cD  is 
assumed  for  analyzing  a  slit  in  a  composite  material.  With  this  assump¬ 
tion  Eq.  (6)  becomes 


Glc  =  A<c  +  Co>Wb  <9> 

Equation  (8)  now  becomes  a  special  case  of  Eq.  (9),  i.e. ,  the  unflawed 
strength  is  recovered  as  c  -*■  0.  This  model  yields  the  following 
relationship  for  residual  strength  o-^  in  terms  of  initial  unflawed  strength 

%  M 


aR  = 


(10) 


5 


t 
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It  is  now  desired  to  derive  an  analogy,  in  terms  of  residual  strength, 
between  the  local  damage  inflicted  by  a  small  hard  particle  impact  and 
damage  inflicted  by  implanting  a  crack  of  known  dimensions  >n  a  static 
tensile  coupon.  When  a  plate  is  subjected  to  a  local  hard  particle  impact, 
the  damage  inflicted  will  be  a  complex  function  of  the  actual  impact  event. 
From  a  practical  standpoint,  however,  it  is  suspected  that  for  impact 
velocities  less  than  the  penetration  velocity,  the  degree  of  damage  is 
strongly  influenced  by  the  amount  of  kinetic  energy  imparted  to  the  plate. 
Such  an  assumption  is  the  basis  of  the  Izod  or  Charpy  impact  test,  and 
is  used  in  developing  the  desired  analogy.  In  particular,  it  is  assumed 
that  the  difference  between  the  energy  density  required  to  break  an 
undamaged  specimen  and  the  energy  density  required  to  break  an  impacted 
specimen  is  directly  proportional  to  the  kinetic  energy  imparted  to  the 
specimen,  Wj^£,  dissipated  over  ome  volume  of  the  specimen.  Thus, 


Ws  *  Wb 


(11) 


where  V  is  an  unknown  volume  over  which  the  kinetic  energy  is  dissipated. 
A  complex  analysis  of  the  exact  impact  event  would  be  necessary  to 
determine  V  theoretically.  It  is  further  assumed,  however,  that  this 
volume  can  be  characterized  by  some  characteristic  surface  area,  A  . 
which  is  independent  of  the  kinetic  energy  of  impact,  and  the  plate  thick¬ 
ness,  t.  Equation  (11)  can  now  be  written  in  the  form 
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Ws-  *b  *  KWKE 


(12) 


or  in  a  more  useful  form 


W, 


=  W  - 
s 


KW 


KE 


(13) 


where 


K*  = 


Wke 

t 


To  find  a  relationship  between  the  impact  damage  and  a  mechanically 
induced  crack,  Eq.  (12)  is  substituted  into  Eq.  (9)  with  the  result 

Glc  r  A(c  +  c0)(Ws-  KWKE)  (14) 

Eq.  (14)  is  now  equated  to  Eq.  (8)  and  the  results  solved  for  c,  yielding 

coK^KE 

c  =  - — —  ns) 

(W8-  kwke) 

Substituting  Eq.  (15)  into  Eq.  (10)  yields  the  desired  relationship  for 
residual  strength  in  terms  of  the  impact  kinetic  energy,  i.e.. 


Equation  (16)  suggests  that  residual  strength  can  be  predicted  as  a 
function  of  kinetic  energy  of  impact  by  executing  two  experiments,  a  static 
tensile  test  on  an  unflawed  specimen  and  a  static  tensile  test  on  a  coupon 
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previously  subjected  to  a  single  point  impact.  If  the  specimen  is 
sufficiently  wide  so  that  the  impact  damage  is  reasonably  well  localized 
the  K  factor  should  be  independent  of  specimen  geometry.  The  value  of 
K  may,  however,  depend  on  laminate  stacking  sequence  and  boundary 
conditions  of  the  experiment  (e.g.,  cantilever  specimen  versus  both  ends 
clamped).  It  should  be  noted  that  if  one  assumes  the  impact  damage  zone 
to  be  analogous  in  residual  strength  to  a  specimen  with  a  c.ircular  hole 
instead  of  a  crack,  the  Bowie  fracture  model  used  by  Waddoups,  Eisenmann, 
and  Kaminski  [6]  will  yield  Eq.  (16).  Thus,  the  form  of  Eq.  (16)  may  not 
depend  on  the  assumed  geometry  used  to  develop  the  analogy.  It  is 
important  to  note  that  the  present  analogy  depends  on  the  validity  of  the 
model  used  to  derive  Eq.  (10).  Data  exists  in  the  literature  [1,6]  which 
indicates  that  such  a  model  has  validity  for  interpreting  composite  laminate 
residual  strength  in  the  presence  of  a  crack. 

Since  one  of  the  key  parameters  in  Eq.  (16)  is  Wg,  theoretical  methods 
which  predict  the  area  under  the  static  stress- strain  curve  are  useful  in 
local  impact  damage  studies.  Such  a  method  for  approximating  Wg  has 
been  developed  by  Petit  and  Waddoups  [7].  This  technique,  including 
example  calculations,  is  discussed  in  the  Appendix. 

For  impact  velocities  equal  to  or  greater  than  the  velocity  for  complete 
penetration  (i.e.,  the  impact  produces  a  clean  hole)  the  residual  strength 
becomes  independent  of  the  impact  event.  Thus,  for  local  impact  damage 
the  residual  strength  never  vanishes,  but  reaches  a  lower  limit  when  the 


8 


AFML-TR-73-309 


impact  produces  a  hole.  Based  on  the  data  in  Figure  1,  this  cut-off 
value  of  <r0  can  be  estimated  by  drilling  a  hole  of  the  same  diameter  as 
the  impacting  particle  in  the  laminate  being  characterized  and  measuring 
the  residual  strength.  In  actual  practice  the  residual  strength  will  drop 
below  this  cut-off  as  initial  penetration,  i.e.,  penetration  without  producing 
a  complete  hole,  will  cause  more  damage  by  removing  material  around 
the  hole  than  will  be  caused  when  the  velocity  is  high  enough  to  produce 
a  clean  hole.  Thus,  there  is  a  range  of  impact  kinetic  energies  for  which 
the  residual  strength  undergoes  a  transition  from  the  damage  model 
represented  by  Equation  (16)  and  complete  penetration  as  represented  by 
a  through  hole.  This  concept  is  illustrated  in  Figure  2. 

For  a  better  physical  understanding  of  the  analogy  represented  by 
Eq.  (16),  the  basic  concepts  are  pictorially  illustrated  in  Figure  3. 
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SECTION  III 

EXPERIMENTAL  PROCEDURE 

Composite  laminates  were  impacted  with  0.  177  inch  and  0.25  inch 
diameter  spherical  steel  projectiles  at  several  velocities  and  then  the 
residual  tensile  strength  of  the  damaged  specimens  measured.  The  speci¬ 
mens  used  in  the  initial  portion  of  this  study  were  6  inch  x  0.  5  inch  straight 
sided  tensile  coupons  with  fiberglass  tabs  bonded  to  the  ends  for  gripping. 
Most  specimens  were  12  plies  thick  (0.  1  inch  -  0.  144  inch  depending  on  the 
material),  and  all  were  [o*,  90 symmetric  laminates.  This  orientation 
was  chosen  for  simplicity,  and  is  representative  of  engineering  laminates. 
The  specimens  were  mounted  as  cantilever  beams  and  impacted  normal  to 
the  surface  at  the  center. 

A  tube  with  a  bore  diameter  considerably  larger  than  the  sphere 
diameter  was  used  to  launch  the  projectiles,  and  a  plastic  sabot  or  cup 
protected  them  from  the  tube  walls  and  provided  a  smooth,  flat  rear  face 
for  efficient  launch  acceleration. 

The  launcher  used  either  compressed  air  or  burning  powder  to  pro¬ 
vide  the  accelerating  pressure.  In  general,  the  compressed  air  assembly 
was  used  to  launch  up  to  velocities  of  about  600  ft/sec  and  powder  was  used 
above  that.  The  launcher  is  capable  of  a  maximum  velocity  of  about  9,000 
ft/sec. 

At  velocities  below  about  500  ft/sec  the  compressed  gases  were  vented 
through  slots  in  the  sides  of  the  launch  tube  and  the  sabot  was  caught  at  the 
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end  of  the  tube  allowing  the  projectile  to  travel  on  alone  without  hindering 
the  experiment.  For  velocities  above  500  ft/sec  the  sabot  often  shattered 
when  abruptly  decelerated,  causing  sabot  fragments  to  follow  the  projectile 
to  the  target  thus  interfering  with  the  experiment.  This  difficulty  was  over¬ 
come  by  slightly  drag  decelerating  the  sabot  to  separate  it  from  the  unim¬ 
peded  projectile  and  then  deflecting  it  away  from  the  flight  path. 

The  method  chosen  to  measure  the  velocity  of  the  projectile  was  very 
simple  but  effective.  The  beams  of  two  inexpensive,  low  power  He-Ne 
lasers  were  directed  across  the  flight  path.  When  the  projectile  interrupted 
the  first  beam  the  high  speed  counter  started.  When  the  second  beam  was 
interrupted,  the  counter  stopped.  The  counter  operating  at  a  known  fre¬ 
quency  gave  the  travel  time  over  a  carefully  measured  distance  between  the 
laser  beams.  A  schematic  of  the  test  set-up  is  illustrated  in  Figure  4. 

High  speed  film  was  used  to  determine  rebound  velocity.  Since  the  pro¬ 
jectile  underwent  very  little  deformation  due  to  impact,  the  kinetic  energy 
transmitted  to  the  specimen  could  be  determined  from  the  rebound  velocity. 

After  the  specimens  were  impacted,  the  damage  was  visually  inspected 
and  then  the  specimens  were  tested  to  failure  in  tension.  A  standard 
10,000  lb  capacity  Instron  was  used  for  tensile  testing  with  a  cross-head 
speed  of  0.05  inches /minute. 
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SECTION  IV 


DISCUSSION 

The  test  procedure  outlined  was  used  to  evaluate  several  composite 

material  systems  to  determine  the  effects  of  constituent  properties  on  the 

impact  response  of  a  composite.  The  fibers  and  resins  evaluated  and  their 

tensile  properties  are  listed  in  Table  1.  Equation  (16)  was  used  to  assess 

the  accuracy  of  the  impact  analogy.  In  particular,  K  was  determined  by 

fitting  the  data  to  the  theory  at  one  value  of  the  kinetic  energy.  A  measured 

value  of  Wg  was  then  used  in  conjunction  with  Equation  (16)  to  predict  values 

of  o^/o'q  as  a  function  of  WKE-  Thus,  an  estimate  of  the  residual  strength 

curve  is  obtained  without  determining  k  or  A  . 

e 

The  experimental  results  indicated  little  or  no  change  in  impact  response 
due  to  variations  of  matrix  properties  for  the  resin  systems  evaluated.  This 
is  illustrated  in  Figure  5  for  Modmor  II  fiber  with  three  different  matrix 
materials  having  various  strains  to  failure.  The  solid  line  represents  theo¬ 
retical  results  obtained  from  Equation  (16)  in  conjunction  with  the  procedure 
previously  described.  Similar  results  were  obtained  with  A-S  fiber  in  both 
epoxy  and  thermoplastic  matrix  materials  (  Table  2).  Although  the 
residual  strength  curves  were  unchanged,  it  should  be  noted  that  the  thermo¬ 
plastic  matrix  composites  displayed  significant  differences  in  failure  mode. 
The  thermoplastic  matrix  composites  displayed  a  denting  mode  of  damage 
with  no  delamination,  while  all  of  the  epoxy  laminates  displayed  local  crush¬ 
ing,  cracking,  and  delamination. 
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Significant  changes  in  impact  response  due  to  variations  of  fiber  proper¬ 
ties  was  apparent.  This  is  clearly  illustrated  in  Figure  6  which  shows  the 
impact  response,  using  Equation  (16)  as  previously  discussed,  of  the  differ¬ 
ent  graphite  fibers  evaluated  and  also  E -glass.  All  of  these  materials  had 
the  same  epoxy  matrix  except  for  the  glass  which  was  Scotchply/ 1002  mater¬ 
ial.  This  data  agrees  well  with  the  results  of  Novak  and  DeCrescente  [8], 
and  Chamis,  Hanson,  and  Serafini  [9].  The  impact  response  of  composites 
appears  to  be  primarily  controlled  by  the  strain  energy  at  failure  of  the 
fibers,  i.e.  the  combination  of  high  strength  and  high  strain  properties  in  a 
fiber  provide  composites  with  higher  resistance  to  impact  damage. 

In  general,  all  of  the  experimental  data  agreed  reasonably  well  with  the 
response  predicted  by  the  model  as  illustrated  in  Figures  7  and  8.  There 
was,  however,  some  concern  that  1/2  inch  wide  specimens  were  too  narrow 
and  that  the  data  was  being  affected  by  edge  conditions,  i.e. ,  the  specimen 
was  not  wide  enough  for  the  damage  area  to  be  classified  as  a  flaw  in  an 
infinite  plate.  Therefore,  a  series  of  A-S/4617  specimens  one  inch  wide 
were  evaluated.  In  addition,  two  laminate  thicknesses  were  evaluated  to 
determine  if  K  was  independent  of  thickness.  A  cursory  examination  of  the 
results  in  Figure  7  shows  that  the  K  factor  is  not  significantly  affected  by 
thickness,  but  is  affected  by  the  smaller  width. 

In  order  to  more  completely  define  a  proper  ratio  of  specimen  width  to 
projectile  diameter,  a  series  of  impact  experiments  were  performed  on 
Scotchply/ 1002  material  using  two  projectile  diameters  and  specimen  widths 
up  to  1-1/2  inches.  Results  in  Figure  8  show  that  a  single  value  of  K  can 
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adequately  describe  the  data  for  W  /d  >5.  It  is  important  to  note  in  both 
Figures  7  and  8  that  the  damage  process  is  slowed  considerably  by  assuring 
that  a  local  impact  occurs.  The  actual  variation  of  K  with  W /d  is  illus¬ 
trated  in  Figure  9. 
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SECTION  V 
CONCLUSIONS 

An  analogy  in  terms  of  residual  strength  has  been  developed  between 
damage  inflicted  by  a  localized  single  point  hard  particle  impact  and  damage 
inflicted  by  inserting  a  flaw  of  known  dimensions  in  a  static  tensile  coupon. 
The  procedure  allows  the  local  impact  resistance  of  laminated  composites  to 
be  characterized  on  at  least  a  comparative  basis  with  a  minimum  of  experi¬ 
mental  data.  The  analogy  is  applicable  to  velocities  which  are  less  than  the 
velocity  of  significant  penetration.  For  penetration  the  residual  strength 
can  be  characterized  as  a  laminate  containing  a  hole  of  the  same  diameter  as 
the  impacting  particle.  Comparison  of  experimental  data  to  the  theory  shows 
good  agreement. 

Several  aspects  of  the  procedure,  however,  need  to  be  pursued  further. 
In  particular,  more  data  is  necessary  to  determine  the  effect  of  impact 
boundary  conditions  on  K.  In  addition,  the  procedure  needs  to  be  applied  to 
laminate  orientations  other  than  [o*,  90°]^  and  to  various  stacking  sequences 
of  the  same  ply  orientations. 
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APPENDIX 

STIFFNESS  AND  FRACTURE  PREDICTIONS  USING 
LAMINATION  THEORY  AND  MAXIMUM  STRAIN  FAILURE  CRITERION 

Lamination  theory  can  be  used  to  predict  the  stiffness  of  a  laminate. 

As  the  laminate  is  loaded  to  failure  individual  ply  failures  occur  until  the 
last  ply  fails.  As  a  ply  fails  the  stiffness  of  the  laminate  decreases  and 
must  be  recalculated  by  deleting  the  ply  which  has  failed  from  the  calcula¬ 
tion.  Thus,  by  using  lamination  theory  to  predict  laminate  stiffness  and  a 
failure  criterion  to  predict  ply  failure,  a  stress  strain  curve  for  the 
laminate  can  be  approximated  and  the  ultimate  strength  and  work  to  break 
predicted.  The  results  of  this  calculation  may  then  be  used  to  estimate  the 
effective  strain  energy  rate  and  fracture  toughness  for  a  Mode  I  failure 
process . 

The  engineering  properties  characteristic  of  an  orthotropic  ply  are  used 
to  determine  the  plane  stress  moduli  of  laminated  plate  theory. 

Engineering  Properties: 

Etj  =  21  x  106  psi 

Ea  =  1 . 7  x  106  psi 

Gj 2  =  0. 6  x  106  psi 

V|  2  —  0  *  28 

V2j  -  V'lZ  £22/^11  =  0.0226 
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Plane  Stress  Moduli: 

Qu  =  Eh  /  ( 1  ■  viz  va  ) 

=  21.2  x  10*  psi 

Qa  =  E 22/(1  -  v  1 2 v 21 ) 

=  1.71  x  106  psi 

Qi2  =  V21Q11  =  0.48  x  106  psi 

Q16  =  0 

Q26  =  0 

Qfct  =  G12  =  0.65  x  106  psi 

The  plane  stress  moduli  for  each  ply  must  be  transformed  to  the 
orientation  of  the  ply  in  the  laminate. 

Transformed  Plane  Stress  Moduli: 

Q11  =  Qu  cos4  0  +  2(Qi2  +  2Q66 )  sin2  0  cos2  0  +  Q22  sin4  0 

Q 2z  =  Qu  sin4  0  +  2(QjZ  +  2Q66)  sin2  0  cos2  0  +  Q22  cos4  0 

Q12  =  (Q11  +  Q22  -  4QM)  sin2  0  cos2  0  +  Q]2  (sin4  0  +  cos4  0) 

Qu  =  (Qii  +  Q22-  2QW)  sin2  0  cos2  0  +  Qji  (sin4  0  +  cos4  0) 

Q16  =  (Qu  _  Q12  "  2Qj$)  sin  0  cos3  0  +  (Qi2  -  Q22 

+  2QW)  sin3  0  cos  0 

<3*  =  (Qu  -  Qia  -  2QW)  sin3  0  cos  0  +  (Qi2  -  Q22 

+  2QW)  sin  0  cos3  0 
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For  the  plies  of  a  (0,  90)g  laminate: 


0 

0° 

90° 

Qn 

2 1 . 2  x 

106 

1.71  x  10 

QZ2 

1.71  x 

io‘ 

21.2  x  10 

Ql2 

0.  48  x 

106 

0.  48  x  10' 

Ql6 

0 

0  ' 

Qab 

0 

0 

Q66 

0 . 65  x 

106 

0 . 65  x  10 

The  stiffness  of  the  laminate  is  obtained  by  summing  the  plane  stress 
moduli  through  the  thickness  in  proportion  to  the  percentage  of  the  thick¬ 
ness  the  k**1  ply  occupies  of  the  n  ply  laminate. 


where 


ak  =  Ahk/h 

tr 

Ah  =  ply  thickness 

h  =  laminate  thickness 
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The  overall  engineering  properties  of  the  laminate  are  obtained  from 
the  inversion  of  the  stiffness  matrix  which  yield  the  laminate  compliances: 


S„ 


S22 


Si  2 


1 

A  22 

E„ 

An  A 22  ■ 

—  2 
A12 

1 

An 

e22 

An  Aj2  - 

—2 

A12 

V12 

^21 

A12 

En 

E22 

An 

—  —2 
A22  "  Aj2 

The  overall  engineering  laminate  properties  for  the  laminate  can  be 
computed  as  the  plies  fail.  As  will  be  shown  in  the  next  section,  the  order 
of  ply  failure  is  90°  and  0°. 

Initial  Laminate: 

—  i  —  0°  l  —  90°  , 

Au  =  2°n  +  I0"  =  II<46  x  10 

-  1-0°  1-90°  ,  . 

A 22  =  2  Q22  +  2  Q22  =  n'46  x  106 

_  i_0°  1-90°  . 

A|2  =  2  °>2  +  z  Q1*  =  °-48  x  10 
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_  i-0°  i-90°  . 

Am  =  1  Qw  +  1  Qu  =  0.65  X  106 

S„  =  S22  =  0.0875  x  10' 6  in2 /lb 

S,2  =  -0.00366  x  10*6  in2 /lb 

Su  =  1.54  x  10'6  in2 /lb 

After  the  90°  ply  fails: 

-  1-0° 

A„  =  ~  Q„  =  10.6  x  106 

1-0 ° 

A22  =  2Q22  =  -86  X  106 

!_0° 

A12  =  -Q1Z  =  .24  x  10* 

—  1  —0°  , 

A-4*  =  -Qw  =  .11  X  106 

which  can  be  converted  to  the  compliances. 

Maximum  Strain  Failure  Criterion 

The  orthotropic  ply  is  characterized  by  six  ultimate  strain  allowables. 
If  any  one  of  the  ultimate  strains  is  exceeded  by  a  ply  of  the  laminate,  the 
ply  has  failed. 
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Ultimate  Strains: 


£  1 

=  0.0085 

-£  l 

=  -0.0098 

£  2 

=  0.0045 

-e  2 

=  -0.0090 

£  6 

=  0.0150 

-£6 

=  -0.0150 

With  the  laminate  under  uniaxial  loading  the  axial  strain,  tx>  which 
causes  failure  in  one  of  the  plies  can  be  computed  by  a  transformation  of 
the  ultimate  strains  for  a  ply. 

«x  =  ci /(cos2  8-  v12  sin2  8} 

£x  =  c2/(sin28-  vJ2  cos2  8) 

cx  =  cj  /(-2  sin  0  cos  6(1  +  vJ2) ) 

The  smallest  axial  strain  which  causes  failure  of  any  one  of  the  plies 
of  the  laminate  determines  the  order  of  ply  failures.  As  each  ply  fails 
the  laminate  stiffness  is  recalculated  to  reflect  the  deletion  of  the  failed 
ply.  The  ply  failure  strains  and  intermediate  laminate  moduli  lead  to  a 
prediction  of  the  stress  strain  curve  and  the  ultimate  strength. 
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Ply  Failure  Strains’  for  0/90  Laminate 

The  90°  ply  fails  first  by  the  positive  transverse  ultimate  strain. 

0  =  90° 

cx  =  £2/(sin20-  vi2  cos2  9)  =  e  2  =  0.0045 

The  0°  ply  fails  last  by  the  positive  longitudinal  ultimate  strain. 
0  =  0° 

£  x  =  £i/(cos20  -  v12  sin2  0)  =  £  i  =  0. 0085 


Laminate  Stress  Strain  Response 


Ply 

Failure 

£x 

E„ 

Aex 

A<rx  =  EAex 

2Ao-x 

0 

o 

O' 

0.0045 

11.46  x  106 

.0045 

51, 525 

51, 525 

0° 

0.0085 

10.6  x  106 

.0040 

42, 000 

93,  525 

.  * .  (tq  =  93,  525  psi 

The  area  under  the  plot  of  ex  versus  ZAo-x  yield  the  work  to  break 
W8  =  455.0  in-lb/in5 

Following  the  arguments  of  Halpin,  Jerina  and  Johnson  [1J,  the 
effective  strain  energy  release  rate  will  be 

Gjc  =  A  c0Wjj  =  169  in-lb/in2 
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for  c  —  0.04  in.  Using  Equation  (1)  wherein  the  S. .  terms  are  the  initial 
o  U 

compliances  for  the  laminate,  the  apparent  fracture  toughness  is 


[Glc1 

1/2 

169  x  10s 

_f(V. 

12.  98 

36.  1  ksiVin 
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Adjusted  for  Rebound  Velocity 
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Figure  1.  Comparison  Between  Strength  Reduction  Produced  B 
A  Drilled  Hole  and  the  Strength  Reduction  Produc 
By  A  Bullet 


Figure  3.  Schematic  of  Residual  Strength  Analogy 
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Schematic  of  Experimental  Set-up 
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Effect  of  Matrix  on  Impact  Residual  Strength 
(MOD  II  Fiber) 


Figure  6.  Effect  of  Fiber  on  Impact  Residual  Strength  (Fiber 
Data  in  Table  1,  ERL  4617  Resin  System  Used  in  All 
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Figure  8.  Effect  of  Specimen  Width  to  Projectile  Diameter 
on  Impact  Residual  Strength  (Scotchply/1002 
Composites) 
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